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Abstract 

Recently, stretched exponential decay of multiple correlations in the periodic 
Lorentz gas has been used to show the convergence of a series of correlations which 
has the physical interpretations as the fourth order Burnett coefficient, a gener- 
alization of the diffusion coefficient. Here the result is extended to include all 
higher order Burnett coefficients, and give a plausible argument that the expansion 
constructed from the Burnett coefficients has a finite radius of convergence. 



1 Introduction 

The Lorentz gas is a model used in statistical mechanics, consisting of a point particle 
moving at constant velocity except for specular collisions with smooth (specifically C 3 ) 
convex fixed scatterers in d > 2 dimensions. The original model || has randomly placed 
scatterers in infinite space, and is thought to have power law decay of correlations, so that 
the Burnett coefficients (denned below as sums of such correlations) are not generally 
expected to exist || [IT| . Here we consider a periodic arrangement of scatterers which 
is equivalent to a dispersing billiard on a torus, for which it is known that two-time 



correlations of the discrete (collision) dynamics decay exponentially j|, [T^]. This, together 
with the finite horizon condition, that is, that the time between collisions is bounded, 
implies the existence of the diffusion coefficient (D^ below). 

A recent paper gives a stretched exponential decay of multiple correlations f5|, and 
uses this to show (again with finite horizon) that the fourth order Burnett coefficient 
(D^ below) exists. Here we extend this result to all the Burnett coefficients. A common 
example for d = 2 with a finite horizon is given by circular scatterers on a triangular 
lattice; for d > 2 the finite horizon condition requires either nonspherical scatterers, or 
more than one scatterer per unit cell. The Lorentz gas and a number of extensions are 
discussed in Ref. fT0|j . 



The Burnett coefficients D^" 1 ' discussed in this paper are defined using series of corre- 
lation functions. Section 2 defines these series and gives three basic results about them. 
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Section 3 gives the main result of this paper, the proof of convergence of these series. 
The series arise in a physical description of diffusion, however the derivation involves 
hydrodynamic approximations and interchange of limits which have not been justified 
rigorously for the Lorentz gas. The physical motivation together with the non-rigorous 
derivation of the series given below from previously stated formulas is given in the final 
section, together with a conjecture about the Burnett expansion. 

The author is grateful for helpful discussion with N. I. Chernov, E. G. D. Cohen, J. 
R. Dorfman and P. Gaspard, and for the support of the Engineering Research Program 
of the Office of Basic Energy Sciences at the US Department of Energy, contract #DE- 
FG02-88-ER13847, and the Nuffield Foundation, grant NAL/00353/G. 

2 Definitions 

In the following, 4>(x) is the billiard map defined on the collision space M, consisting 
of points x = (r, v) G M for which the position r G R d is on the boundary of one 
of the scatterers and the velocity following a collision v G R d is of unit magnitude 
in an outward direction from the scatterer. Greek indices a,/3, ... = l,...,d denote 
components of vectors and tensors in R d , and a dot a • b denotes the usual inner product 
J2 a a ab a corresponding to the Euclidean metric. We have two functions T : M — > R and 
a : M — > R d which describe the embedding of the collision dynamics into physical space 
and time, as follows. T(x) is the time (also distance since the speed of the particle is one) 
between the collision at x and the next; it is a piecewise Holder continuous function fl], ||. 
&(x) is the lattice translation vector associated with this free flight when the configuration 
variable r is unfolded onto a periodic tiling of R d ; it is a linear combination of the lattice 
basis vectors with integer coefficients, and is a piecewise constant function. The 
finite horizon condition ensures that both T and a are bounded. The average (■) denotes 
integration over M with respect to the invariant equilibrium measure. In terms of this 
average we define AT : M -> R by AT(x) = T(x) - (T) so that (AT) = 0. 

The billiard dynamics is time reversal invariant, that is, there exists an involution 
T : M — > M (given simply by the specular reflection law) with the property 



<fi o T o (p = T 



(1) 



In addition, T preserves the equilibrium measure, that is, 



(Tog) = (g) 



(2) 



for arbitrary measurable function g : M — > M. The map T also satisfies 



T o T o (p 
a o T o (p 



T 



—a 



(3) 
(4) 



Thus (a) = 0. 
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The wave vector k is to be understood as a formal real expansion parameter with d 
components (although physically we would like to interpret it as a vector with a value in 
R d ). The dispersion relation s[k] is to be understood as a formal power series 

oo 

»M= E 1 " 1 E D^l a j ai ...k am (5) 

m=2 ot\...a m 

in terms of the Burnett coefficients which are assumed to be real, totally symmetric 
tensors of rank m. That is, an equation (specifically Eq. (|lf|) below) involving s[k] is 
to be interpreted as a sequence of equations (specifically Eq. (|16|) below) obtained by 
equating coefficients of powers of k. The symbol 1 denotes 

The existence of Burnett coefficients satisfying the equations flUf ) is not assumed a 
priori; we show in Lemma [l] below that equations fll~6|) express the 1)/2 independent 
components of .D^ as series not containing any of the D^ m \ then the d(d + + 2)/6 
independent components of as series containing only the and so on. Lemma |2] 
shows that they are indeed real, and Thm. |] shows that the limit exists. 

We define formal power series / and F by 

f[k] = s[k]AT + ik-a (6) 

n-l 

F i k ] = E /W ^ (7) 

i=— n 

where the dependence on x and on the positive integer n is suppressed in the notation; 
the limit n — > oo will be taken later. We have (/) = and (F) = at each order in k 
and for each n as a consequence of (AT) = and (a) = above. 

We define cumulants Qn^\ (also formal power series) for integers N > 2 as 

Qiv[kj- ^ (-1) fjT^TTKl (8) 

with j and ^ integers satisfying j > 2 and i/j > 0, and z/ = X)j v j is the total number of 
correlations in the product. For example 

Q2 = (F 2 )/2 (9) 

Qs = (F 3 )/6 (10) 

Q 4 = ((F 4 )-3(F 2 ) 2 )/24 (11) 

g 5 = ((F 5 )-10(F 3 )(F 2 ))/120 (12) 

Q 6 = ((F 6 )-15(F 4 )(F 2 )-10(F 3 ) 2 + 30(F 2 ) 3 )/720 (13) 

Now Qn contains exactly N powers of F, and so it contains terms k m only for m > N, 
and we can write it as 

oo 

Q N [k] = E Qn (14) 

m=N a.\...a m 
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thus defining totally symmetric tensors qN, m for m > N. 

The Burnett coefficients are found by equating the formal power series on both sides 

of 

1 oo 

s[k] = i^^yE^[k] (15) 

that is, 

i m 

mjj(m) = n y / lg N 

ai...a m n— -»oo 2n(T) ' 

These equations determine the explicitly as real tensors, subject to convergence of 
the limit, as shown by the following two lemmas. 

Lemma 1 The right hand side of Eq. ftldj ) does not contain Z^" 1 ') such that ml > m. 

Proof: We have N > 2, and each Qn contains N powers of F, thus each term has at 
least 2 powers of F. Each F has at least 1 power of k, and there are m powers of k in 
total, so each F has at most m — 1 powers of k. D^ m ) appear in F associated with ml 
powers of k, so ml < m — 1 for any Z)( m ') appearing. 

Remark: It is possible there are no factors of D^ m ) on the right hand side, in fact the 
lemma shows that this is true for m = 2. The case m = 2 can easily be written explicitly; 
Eq. flTBD becomes 

i n— 1 n— 1 



This is a discrete time version of the well-known Green-Kubo formula for the diffusion 

(2) 

tensor (which reduces to a single diffusion coefficient in the isotropic case D a p = D5 a p). 
An equivalent discrete time equation appears in , also for m — 4. 

Lemma 2 Despite the appearance of the imaginary number i in the above definitions, 
the Burnett coefficients are real if they exist. 

Proof: We note from the definitions that s[k], /[k] and F[k] have pure imaginary 
coefficients for odd powers of k and real coefficients for even powers of k. This property 
is preserved by addition and multiplication of power series, so it also holds for the Qjv[k]. 
This implies that the qN, m are imaginary for odd m and real for even m. The result 
follows from Eq. (|TB|). 

Before proceeding with the more technical convergence proof, we note another im- 
portant result: 

Lemma 3 _D( m ) = for m odd. 

Proof: From the properties of the time reversal operator T given above, (F J ) has 
zero contribution from any term with an odd number of a factors. The result follows by 
induction on m: assume that D^ m ' = for all odd ml < m, then by Lemma [I] all terms 
in s[k] contributing to Z)( m ) have even powers of k, and from the oddness of a under time 
reversal, so also do the lk • a terms. Thus D^ m \ which is constructed from terms with m 
powers of k, must be zero for m odd. 
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3 Convergence of the series 



The averages (F^) appearing in the cumulants contain summations over j variables with 
range — n to n — 1, and could grow as fast as 0{tiP) in general. Thus each term, which 
is a product of such averages could grow as 0(n N ) in general. For the limit in Eq. flT6|) 
to exist, we require that the series grows only as 0(n). Although the growth of each 
product of correlations cannot be controlled this well, cancellations occur in constructing 
the cumulants. This is expressed in the following theorem which, together with Lemmas [l] 
and 0, implies the existence of the Burnett coefficients: 

Theorem 4 qN, m is defined in Eqs. (fj|, [|, \L\ ) for integers N and m satisfying 2 < N < 
m. The limit 

lim -qN,m;ax...a m (18) 
n — >oo IX 

exists for all such N and m in the periodic Lorentz gas. 

The structure of the proof of Thm. |] is as follows. We state the theorem expressing 
stretched exponential decay of multiple correlation functions. Next, the terms appearing 
in ( |l6l) are written as a time ordered sum, so that this theorem can be applied. Then 
we show that all the terms connected by the application of the theorem have coefficients 
which sum to zero, so that only the stretched exponential corrections remain. Finally, a 
bound of n multiplied by a polynomial is put on the number of terms at each order of 
the stretched exponential, so that the series divided by n converges absolutely. 

Thm. |] is based on the following result: 

Theorem 5 (Theorem 2 of Ref. j^J) Let %\ < ■ ■ ■ < and 1 < t < k — 1. Then 

Kfi 1 ■ ■ ■ fi k ) - (fi 1 ■ ■ ■ ft)(ft+ + i ■ ■ ■ fi k )\ < c k ■ \h - il \ 2 \^-^ l/2 (19) 

where Ck > depends on the functions fx, . . . , fk, and X < 1 is independent of k and 

fli ■ ■ ■ i fk- 

The theorem applies to piecewise Holder continuous functions fj such that (fj) = for 
all j and uses notation /j = fj o0 ! . As noted in Ref. ||, we expect based on Refs. [0, f| 



that it should be possible to prove a stronger bound A'* t+1 but the above bound is 
sufficient for our purposes here. 

The qw m as defined in the previous section are finite sums of terms of the form (see 
Eqs. (0J, 0)) 



£ ,; l) " 'ifoTw £ ->.-<- ^> (20) 



i\...i N - 



multiplied by constants such as the lower order Burnett coefficients. The / here and for 
the remainder of this section are T or a, both of which satisfy the conditions of Thm. [5[ 
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The exact number of terms of this kind is not important; it depends on N and m but 
not n and therefore does not affect convergence of the limit n — > oo. 

In order to use Thm. [5] we need to put the times i p in numerical order. The unre- 
stricted sum over all the i p is replaced by an ordered sum i\ < z 2 ...ijv over all N\/S[i] 
permutations of the i p . S[i] is a symmetry factor to account for the fact that some of the 
i p may be equal; the exact form is unimportant since it is a common prefactor, indepen- 
dent of the Uj. Not all N\ permutations of the correlations are distinct: it does not matter 
in which order the fj are multiplied within a correlation, or which order correlations of 
equal numbers of fj are multiplied; thus both factorials in the denominator disappear, 
leading to 



E 



(-irv-i)! 

+ permutations}] 



il<i2...i N 



(21) 



The "permutations" remaining in (pTJ) consist of the remaining N\/([\j Ujljl"*) — 1 rear- 
rangements of the i p that are not equivalent by reordering the product of correlations or 
the product of / within a correlation. 

As an example, we give the expression for N = 6: 



E 



J / / Hi fi2 Hz Ha f is He \ 

-i r-n \\Jl J2 J3 J i J5 J6 I 



S\i) 

- \ifi muv. iv. i -rm + (ft ft) (ft ft & ft) + uvmiJ7JVJ-:f«) 
■iii iv) (J7. iv. iv m + (ftft)(ftftftft) + (ft ft) (ft ft ft ft) 
+ (n 2 ft) (ft ft n 5 ft) + (ftft)(ftftftft) + (ftft)iftftftft) 
+(ftft)(ftftftft) + (ft ft) (ft ft ft ft) + (ftft)(ftftftft) 
+(ftft)(ftftftft) + (ftft)(ftftftft) + (ftft)(ftftftft)\ 

- '(ft ft ft) (ft ft ft) + (ft ft ft) (ft ft ft) + (ft ft ft) (ft ft ft) 

+(ftftft)(ftftft) + (ft ft. 17) (.17. IV ft) + (ft. 17. 17) (.17. 17 ft) 
7.17. 17. 17) (.17. 17. 17) + (ft ft 17) (ft 17 17) + (.17 .17. 17) (.17. 17. 17) 
■ (ft. 17. 17) (.17. 17. 17) 

+2 



(22) 



(ftft)(ftft)(ftft) + (ftft)(ftft)(ftft) + (ft ft) (ft ft) (ft ft) 

7.17. 17) (.17. 17) (.17. 17) + (.17 .17) (.17. 17) (.17. 17) + (ftft)(ftft)(ftft) 
■(ft ft) (ft ft) (ft 17) + (.17. 17) (.17. 17) (.17. 17) + (ft ft) (ft ft) (ft 17) 
7. 17. 17) (.17. 17) (.17. 17) + (.17 .17) (.17. 17) (.17. 17) + (ftft)(ftft)(ftft) 
7.17. 17) (.17. 17) (.17. 17) + (ft .17) (.17. 17) (.17. 17) + (ft ft) (ft ft) (ft ft)}} 

Here, the four terms correspond to the partitions of 6 which do not contain 1; in the 
above notation the nonzero Uj are {uq = 1} with 6!/6! = 1 term; {i>2 = 1, z/4 = 1} with 
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6!/2!4! = 15 terms; {u 3 = 2} with 6!/2!3! 2 = 10 terms; and {v 2 = 3} with 6!/3!2! 3 = 15 
terms; compare with Eq. fll3|) . 

Now we apply Thm. [5] to the largest gap, i t +i — it- Any of the largest gaps will suffice 
if more than one is largest. Before tackling the general case, we see how it works in 
the N = 6 example. Notice that, whatever the value of t, the theorem combines all the 
above correlations to leave terms (the number of which is a function of N) bounded by 
X' H+1 ~ lt ' multiplied by powers of the time differences. Explicitly, for t = 1, all terms 
cancel individually because (fj) = 0. For t = 2 the (/ 6 ) term cancels with one of the 
(f 2 )(f 4 ) terms, six other (f 2 )(f ) terms cancel with three of the (/ 2 ) 3 terms and the 
remaining terms all split leaving an (/) term. For t = 3 the (f 6 ) term cancels with one 
of the (/ 3 ) 2 terms, and all of the others split leaving an (/) term, t = 4 is analogous to 
t — 2 and t = 5 is analogous to t — 1. 

In general we must show that the coefficient (— l)" _1 (z/ — 1)! in Eq. ([H]) combined 
with the numbers of terms of various types leads to complete cancellation for all values 
of N. Consider a general term (ignoring the S[i] which is the same for each term) which 
is unaffected by a split at time t. Each correlation contains times i p < t or times i p > t 
but not both. Thus it can be written schematically as 



where all times i p to the left of the bar "|" are less than or equal to t and all times to 
the right of the bar are greater than t. Let there be A correlations to the left and B 
correlations to the right, so A + B = v. 

This term will cancel (up to stretched exponential corrections) with any term which 
is split to the same form, if the sum of the coefficients (the (— \y~ x (v — 1)!) is zero. The 
terms that are split to a given form consist of correlations that are either the same as the 
above, or are joined in a pairwise fashion with a correlation on the other side of the bar. 

Again, an example is helpful: When iV = 8, a split at t — 4 combines the following 



terms: -GiffffjiftftjlififtHftft) with 2(/f f 2 2 f 5 5 f?)(ft f 8 8 ), 

-(ftffftf^iftftffft) and -(ftf?fi T ft)(ftftftft)- These all cancel because 
-6 + 2 + 2 + 2 + 2-1-1 = 0. 



The term given in Eq. (|23|) has coefficient (— l)^ _1 (z/ — 1)!. There are AB terms with 
coefficient (— l)^~ 2 (z/ — 2)! obtained by combining a single correlation on the left and the 
right. There are A(A - 1)B(B - l)/2! terms with coefficient (-l)^ 3 (z/ - 3)! obtained 
by combining two correlations on the left and the right, and so on until all min(A, B) 
correlations on the side with the fewest correlations have been combined. The total 
coefficient is thus given by 



()()••• 0100 •••() 



(23) 




(24) 



To show that the coefficients cancel, we therefore need the following lemma: 
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Lemma 6 H(A, B) = for all positive integers A and B. 

Proof: The sum is symmetric in A and B so suppose that A > B without loss of 
generality. Then the summand is the product of a constant (— 1) A+B ^ 1 A\, an alternating 
binomial of degree B, that is, (— 1)~ P B\/((B — p)\p\) and a polynomial in p of degree 
B — 1, that is, (A + B — p — 1)!/(A — p)\. We will use summation by parts to lower the 
degree of both until the result is zero. 

We note the summation by parts formula 

B B B B 

^2x p y p = y J2 x P + J2(y<i-y<i-^J2 x p ( 25 ) 

p=0 p=0 q=l p=q 

which can be demonstrated by collecting terms on the right hand side. Now substituting 
— l)~ p £> !/((£> — p)\p\ and y p = (A + B —p — 1)1/ (A — p)\ we can show by induction 



p 




on q from B downwards that 

L ) (B- q y.(q-iy. 9 > u ( 2 q) 
q=0 v ; 

hence the first term on the right hand side of Eq. (KH) vanishes. We can also simplify 



so Eq. now reads 



JA + B-q-l)\ 



H(A, B) = - B) ± <- 1 )'Y - % (A± — - « - (28) 

Shifting the summation index by one we find 

H (A, B) = (1 - B)H(A, B - 1) (29) 

The proof of Lemma ^| follows by noting that H(A, 1) = 0. 

We now conclude the proof of Thm. [|. Recall that the series (|20|) have been rewritten 
in the form (pip. Thm. ^| is applied to (one of) the largest gap(s) Ai max = i t+ i — i t , 
partitioning the terms into subsets which split into a particular form (p3|). Lemma ^| 
shows that the coefficients of all terms in a subset conspire to cancel, so that each subset 
is bounded by the error term in theorem |5], that is A' Almax ' 1/2 multiplied by a polynomial 
in the time differences. 

Finally we estimate the number of terms with each value of Ai max . The first time i\ 
varies freely from — n to n — 1, having a total of In values. One of the time differences 
is equal to Ai max , and the other k — 2 time differences can range from to Ai max , so the 
total number of terms with a given Ai max is less than 2n(k — l)Ai^~^\ in particular a 
polynomial in Ai max multiplied by n. Thus the series divided by n appearing in Thm. ^ 
is bounded by a product of polynomial factors and the decaying stretched exponential, 
and hence converges absolutely. This concludes the proof of Thm. (| and the proof of 
existence of Burnett coefficients. 
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4 Physical motivation and remarks 



This section makes the connection between the Burnett coefficients defined in the pre- 
vious sections and equations found in the physics literature. The latter equations are 
phenomenological and have not been shown rigorously in a limiting fashion from the 
Lorentz gas, and a few nonrigorous limit interchanges are made to connect them with 
the expressions defined in the previous sections. First we consider the dispersion relation, 
the equation for the Burnett coefficients, and finally the whether the dispersion relation 
can be used to define an analytic function. 

The dispersion relation (|5|) with k interpreted as a real vector represents the solution 
of a generalised diffusion equation proposed by Burnett || containing higher derivative 
terms that become important on small scales, 

oo 

dtP= E E D^ am d ai ...d amP (30) 

m=2 ai...a m 

assuming a solution of the form 

p(r,t) ~ exp(s(k)t + lk ■ r) (31) 

Here, d a = d/dr a . Nonlinear terms such as powers of d a p are excluded on physical 
grounds since p is a projection onto real space (r e R d ) of a phase space density satisfying 
a linear evolution equation. The phase space is a subset of R 2dM corresponding to the 
possible positions and velocities of M > 1 particles. The dispersion relation is a more 
robust formulation than the generalized diffusion equation (^) since the former may be 
supplemented by nonanalytic functions of k to account for situations (other than the 
periodic Lorentz gas) in which some of the Burnett coefficients do not exist. 

Chapter 7 of Ref. [|7| obtains the dispersion relation from the microscopic dynamics 
using the equation (7.91 in this reference): 

n-l 

1 = fim ( J] exp[-s(k)T(^x) - ik • a(0 i x)]) (32) 

i=—n 

We write T = (T) + AT as in previous sections, take out the constant factor of (T), and 
take the logarithm to find 

Kk)= ffin o ^yMexp[F(k)]) (33) 

where F is defined (as a power series) in Eq. (|7|). Now the exponential and the logarithm 
are expanded in power series and the resulting terms containing iV powers of F are 
collected to become the cumulants Qn defined in Eq. (||). The cumulant form of the 
expansion is possibly more robust than the above equations due to the cancellations 
among the terms that combine to construct each cumulant. 

Since it is desirable from a physical point of view to interpret k as a real variable, we 
conclude with the following conjecture: 
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Conjecture 7 The series (pi) converges when ksDc R d /or some nontrivial domain 
T>, and so defines a function s(k) in this domain. 

Note that s(k) (if it exists) is a real function as a consequence of Lemma [^, and 
physically is expected to be negative except at the origin (otherwise the density p would 
grow exponentially with time); this puts further constraints on the Burnett coefficients. 

Unfortunately the proof given in the previous section contains many undetermined 
functions of k, and the Burnett coefficients are defined by a complicated recursive rela- 
tion (0), so a proof is unlikely using the techniques of this paper. 

There are two results that make such a result plausible. The first is that in the 
Boltzmann limit of a hard sphere gas, that is, a gas with many moving particles at 
low density and with recollisions ignored, the expansion in k (in this context called 
the linearized Chapman- Enskog expansion) converges ||. Of course, the hard sphere 
collisions are similar to that of the Lorentz gas, but recollisions cannot be ignored in 
general. 

The second result is exact, but for a highly simplified (piecewise linear) system. We 
consider the map <f> : 1R — > R given by 

<t>{x) = ^ -2x + 3[x] (34) 

where [x] is the greatest integer less than or equal to x. The dynamics defined by <f> is 
equivalent to a random walk where the particle moves with equal probability from one 
interval I„ = (n — 1/2, n + 1/2) to the left, J n _i or to the right, I n +i- The dispersion 



relation s(k) follows directly from the above phenomenological solution (31), 

p(n, t) = exp(s£ + ikn) (35) 

After one iteration, 

p(n, 1) = -[exp(ik(n — 1)) + exp(ifc(n + 1))] (36) 
= cos kexp(ikn) (37) 

leading to 

s(k) = In cos k (38) 
which has a power series around k = with a radius of convergence equal to ir/2. 
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